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We show that renormalization group (RG) theory applied to complex networks are 
useful to classify network topologies into universality classes in the space of configura- 
tions. The RG flow readily identifies a small- world/fractal transition by finding (i) a 
trivial stable fixed point of a complete graph, (ii) a non-trivial point of a pure fractal 
topology that is stable or unstable according to the amount of long-range links in the 
network, and (iii) another stable point of a fractal with short-cuts that exists exactly 
at the small- world/fractal transition. As a collateral, the RG technique explains the 
coexistence of the seemingly contradicting fractal and small-world phases and allows to 
extract information on the distribution of short-cuts in real-world networks, a problem 
of importance for information flow in the system. 
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A generic property that is usually inherent in scale- 
free networks but applies equally well to other types of 
networks, such as in Erdos-Renyi random graphs, is the 
small-world feature [HQ. In small- world networks a very 
small number of steps is required to reach a given node 
starting from any other node. This is expressed by the 
slow (logarithmic) increase of the average diameter of 
the network, f, with the total number of nodes No, f ~ 
IniVo, where r is the shortest distance between two nodes 
through network links. 

The small- world property has been shown to apply in 
many empirical studies of diverse systems. However, re- 
cent work [S-Q showed that many networks that have 
been found to display the small-world property, such 
as the WWW, are indeed fractal, indicating a power- 
law dependence of the distances with the network size, 
f ~ Nl ,dB , where (is is the fractal dimension, up to a 
certain length-scale before the global small- world behav- 
ior is observed. Therefore, it is not clear how it is possi- 
ble that fractal scale-free networks coexist with the small 
world property. This shows the need for a mathematical 
framework that reconciles these two seemingly contradic- 
tory aspects, fractality and the small-world property. 

In this paper we show that renormalization group (RG) 
theory, previously developed to understand critical phase 
transitions in physical systems |7[ and recently extended 
to inhomogeneous networks [3- 6] , provides such a frame- 
work. The main result of our work is four-fold: (1) We 
introduce a method based on the RG technique that clas- 
sifies network topologies into three universality classes ac- 
cording to fixed points of the RG flow. We find a stable 
trivial fixed point of a complete graph, a non-trivial fixed 
point of a fractal structure that becomes stable or unsta- 
ble according to the amount of long-range links added 
to the network, and a third stable fixed point that ex- 
ists exactly at the small- world/fractal transition consist- 
ing of a fractal with short-cuts. (2) The RG technique al- 
lows for finding the distribution of short-cuts overlaying 
a pure fractal topology, a technique that we test in real- 
world networks like the WWW and biological networks. 
(3) The RG identifies a second point which is associated 



with information flow in the system. (4) The RG anal- 
ysis finds an explanation for the seemingly contradiction 
between the small- world effect observed at a global scale 
in real- world networks and the fractal behavior occurring 
at finite scales. 
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FIG. 1: (a) Top panel: fractal network Go with box cover- 
ing for b = 2, and its renormalized network Rt(Go). Bottom 
panel: short-cuts are added (black links in lower left panel) 
we obtain network G'. We apply the RG transformation to G' 
and obtain Gb (lower right panel), (b) Sketch of the number 
of boxes versus the box diameter according to the probability, 
P(r) — Ar~ a , to add short-cuts between nodes at distance r. 
When the network is in the stable phase, s = a/ds > 2, RG 
flows toward the pure fractal fixed point Go, where the num- 
ber of boxes is a power-law for all values of b (red straight 
line). When s < 2, the number of boxes displays a power- 
law with an exponential cut-off at large values, indicating 
that globally the network topology behaves as a small-world, 
but for small values of b the network still exhibits fractality. 
For small values of a, the number of boxes decays exponen- 
tially and the network displays the small-world property at 
all scales. 



We apply the RG to complex networks using the box- 
covering technique [|| . The network is covered with boxes 
such that all nodes within a box are at a distance smaller 
than b (top panel of Fig. UK), where distance is the num- 
ber of links along the shortest path between two nodes. 
Once the network is tiled, we construct the renormal- 
ized network: by replacing each box with supernodes (or 
renormalized nodes) and these supernodes are connected 
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if there is at least one link between two nodes in their 
corresponding boxes. When this RG transformation, R b , 
is applied to a network Go, it leads to a new network Gb- 
If Go is self-similar, as it was empirically shown to be the 
case for the WWW and many biological networks 0-Q , 
the RG leads to a structure that presents similar proper- 
ties as Go- More technically, if Go is a fractal network, 
then i?b(Go) = Go and Go is a fixed point of the RG 
flow. 

Suppose we start with the fractal network Go and add 
short-cuts according to the distance r between nodes with 
probability p(r) = Ar~~ a , where r > 1. The new net- 
work with short-cuts, G', is not self-similar anymore or 
in other words, Rb(G') ^ G' (see Fig. QJi for a simple 
example of this process). Here we show that depending 
on the value of the exponent of the short-cuts a, the ap- 
plication of the RG process brings G' either back to the 
original self-similar structure Go or transforms it into a 
complete graph (where all nodes are connected to each 
other). Go and the complete graph are both fixed points 
in the space of networks with a substantial difference be- 
tween them. Go is an unstable fixed point of Rb since 
a small perturbation (small number of short-cuts) may 
lead it to a drastically different network under Rb- The 
complete graph is a stable or trivial fixed point because 
any small perturbation always returns the network into 
the complete graph under R b . 

We start by calculating the RG flow in the space of con- 
figurations. Let oIb be the box (or fractal) dimension of 
the original self-similar network Go- Thus b dB — No/N b 
is the average number of nodes in a box of size b where N b 
is the number of nodes in Gb (or number of boxes in Go) 
and No the number of nodes is Go- After a rcnormaliza- 
tion step is applied to the network G', the probability to 
find a short-cut between two nodes at distance r in the 
renormalized network Gb (black links in lower right panel 
of Fig. QJi) is Pb(r) = 1 — (1 — p{br)) h B , and therefore, 



Pb (r) = 1 - [1 - A{brY 



(1) 



For simplicity we write x = A (br) a and B(r) = 
j^2d B /a r -id B _ After repeatedly applying the renormal- 
ization transformation, i.e. b — > oo, we find a fixed point 
of the RG flow defined at 



p*(r) = lim pb(r) = 1 - 

b— >oo 
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(2) 



Analysis of Eq. ([2]) reveals the existence of a critical 
value ats = a/ds=2 separating two phases of the 
RG flow. If s > 2, we find p*(r) = 0. Therefore, the 
RG flow brings the network toward the self-similar fixed 
point Go, implying that the added short-cuts disappear 
under the renormalization flow. On the other hand, if 
s < 2, we find p*(r) = 1. In this case, G 1 flows under 
RG toward a trivial fixed point consisting of a complete 



graph where all the nodes are connected to each other. 
If s = 2, G 1 flows toward another non- trivial stable fixed 
point consisting of the original fractal network Go with 
short-cuts following p* (r) = 1 - cxp(— Ar~ 2ds ) (Fig[2k). 
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FIG. 2: (a) RG flow diagram, (b) Phase diagram of the three 
phases found via RG at s < 1, l<s<2, and s > 2. 

To better understand the features of the phases identi- 
fied by the RG flow, we analyze the behavior of the aver- 
age network degree under renormalization. This calcula- 
tion allows one to identify a second critical point within 
the unstable phase related to information flow in the sys- 
tem. Let zq be the average degree (number of links per 
node) of the unperturbed network, Go, and z' the aver- 
age degree of G' after the short-cuts are added. Then, 
z' — zq = 2 No ) wnere M(L) is the number of short-cuts 
at distance L (the diameter of Go). Since Go is fractal, 
we find, 

M(L) « d B [ Ar^r^^dr = -^-(L^C 1 -*) - 1). 
Ji 1 — s 

(3) 
(4) 



Hence we obtain 
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After renormalizing the network G 1 with length-scale b, 
short-cuts connecting nodes inside a box will not appear 
in the renormalized network, Gb- Therefore, the number 
of remaining short-cuts in G b is simply the number of 
short-cuts that connect different boxes, i.e. M(L)—M(b). 
If Zb is the average degree of the renormalized network, 
Gb, then 

2(M(L) - M(b)) , 
z h - z = = (z - zo)f N (b), (5) 



where 



f N (b) = 



N b 



j^d B (l-s) _ fod B (l- 
L<1b(1-s) _ 1 




(6) 



We define the renormalization parameter Xb = N^/Nb = 
b ds , and in the limit of large networks, L — > oo, we find 
the scaling: 



fN(x b ) 



• b • 



(7) 
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where the RG exponent A depends on the long-range ex- 
ponent a as 

Equation ([SJ (see Fig. ^jp) identifies two transitions 
separating different phases in the space of configurations, 
as depicted in the phase diagram of Fig. [2] The first tran- 
sition at s = 2 corresponds to the point when A = 0, and 
separates a stable phase with A < for s > 2 from an 
unstable phase with A > for s < 2. Therefore this tran- 
sition corresponds to the complete graph/fractal transi- 
tion identified by the analysis of Eq. ((5} , and corresponds 
to the point at which the network topology dramatically 
changes. In the unstable phase, s < 2, the average degree 
increases (A > 0), so that under infinite steps of the RG 
procedure the network becomes a complete graph with 
infinite average degree in the thermodynamic limit. On 
the other hand, when s > 2, the the network conserves 
the global fractal structure of Go . Under the RG flow the 
difference between Zb and zq goes to and the short-cuts 
disappear, returning G back to its original fractal struc- 
ture. In this state, the diameter of the network grows 
as a power law with the number of nodes, implying a 
large-world fractal structure. Instead, if long-range con- 
nections are added with s < 2, the small- world property 
is achieved, where the diameter of the network grows log- 
arithmically with the number of nodes. Therefore, the 
s = 2 (or a = 2ds) transition is a small- world/fractal 
transition. This calculation generalizes the small/large- 
world transition, previously found in Refs. [8l-[To| for lat- 
tices to the case of complex networks. 

An important point readily emerges from the analysis 
of Eq. (|5]) when s = 1 (a = ds) within the unstable 
phase, as shown in Fig. [2^. Notice that s = 1 coincides 
with the optimal point of decentralized navigability of 
Kleinberg [ll| for lattices with fractal dimension ds |12l |. 
and therefore this results could be seen an a plausible 
extension of the results in Refs. [IJ Ell f° r scale-free 
complex networks. We expand on this point later. 

As a test of the RG predictions we use a model of frac- 
tal networks, as described in Ref. Using the fractal 
model, short-cuts with an exponent a can be added to 
the network and the prediction of Eq. ([8]) can be tested 
in a controlled manner. 

The fractal model network is built as follows [5\: At 
generation n = 0, we start with a star network of 5 nodes, 
i.e. a node in the center and four nodes connected to the 
center node. Then, generation n + 1 is obtained recur- 
sively by attaching to new nodes to the endpoints of each 
link I of generation n. In addition, we remove links I of 
generation n and add x new links connecting pairs of new 
nodes attached to the endpoints of I (see top panel of 
Fig. Eh for an example at with n = 2, to — 2, and x = 1). 
The algorithm leads to a pure fractal scale-free network 
with degree distribution exponent 7 = l+ln(2m+a;)/ln 2 
and fractal dimension ds = ln(2m + x)/ln 3. 




FIG. 3: Renormalization applied to the fractal model, (a) 
Average degree of the renormalized network versus Xb for the 
fractal model with n = 6, m = 2,x = 1. (b) Values of A 
obtained from Fig. [3^ (green data points) and the prediction 
of the RG from Eq. (J7J) (solid line), (c) Finite size effects 
in Zb — zq. (d) Navigation time T(a) for the fractal network 
model. Circles show the minimum value. 



Figure [3^ shows the results of the RG flow applied to 
the fractal model network with n = 6, to = 2, x = 1 with 
ds = 1.46 for various long range exponents a. Starting 
from a pure fractal topology we pick a node and add a 
random connection to another node at distance r accord- 
ing to the probability p(r) = Ar~ a when r > 1 (in our 
results we repeat this process for 10% of the nodes in 
the network, although our conclusions are independent 
of this value.) The renormalization is performed numer- 
ically using the box covering algorithm called MEMB in 
Ref. [H. Notice that the MEMB algorithm leads to net- 
works that are smaller than the original, and therefore 
few points are obtained in the plot of 2% — zq vs. Xb- 
To overcome this problem and obtain better resolution, 
we take advantage of the self-similar aspect of the net- 
work and perform a "partial renormalization" , in which 
parts of the network are subsequently renormalized into 
supernodes. We follow the behavior of 2% in the RG flow 
for a given a. For s < 1, the average degree follows a 
power-law with exponent A = 1 as in Eq. ©. When 
s > 1 the exponent follows the theoretical prediction 
Eq. ((SJ, A = 2 — s. Figure [3p shows a very close compar- 
ison with theory indicating the transition between the 
stable region and the unstable region and the optimal 
navigability point. In Fig. |3t we show the dependence 
of Zb — zq for different values of No. When b is large, 
finite size effects become evident and the average degree 
of the renormalized network deviates from the expected 
power-law of Eq. ([7]). 

The fractal model also allows us to verify directly that 
the point s = 1 could be regarded as a plausible extension 
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of the optimal point of decentralized navigability found 
by Kleinberg in lattices [ll|. We test this by measuring 
numerically the average time, T(a), for a message to be 
delivered from a source node to a target node along the 
links of the network. It is important to notice that since 
scale-free networks are not embedded in any euclidean 
space, one cannot directly apply the decentralized algo- 
rithm as introduced by Kleinberg. In the case of scale- 
free networks, we allow for the message holder to have 
information on the distance between any node and the 
target in the fractal background Go, but not on their 
long-range short-cuts that exists in G'. In Fig. [3J: we 
show simulation results for T(a) versus a for the fractal 
model with m = 2, x — 1 and for different values of sys- 
tem size Nq. We find that the value a c (No) correspond- 
ing to the minimum delivery time for a given No slowly 
converges to the critical value a c (No) — > a c = ds = 1-46 
as iVo — !• oo, implying a navigability transition at s = 1 as 
predicted by the RG analysis. A finite size scaling anal- 
ysis of a c (N ) versus the inverse square of the logarithm 
of N , as suggested in Ref. [HI, confirms this result. 

An advantage of the RG approach is that it allows 
for a measurement of the type of short-cuts present in 
real-world networks. Previously in this paper, we started 
with a pure fractal structure Go to which short-cuts were 
added, generating network G', and analyze the stability 
of G' under the RG procedure. We now switch to the 
study of real- world networks where we tackle the inverse 
problem. The real- world networks we examine are known 
to have an underlying fractal structure since the mea- 
surement of Nb versus b leads to a power-law relation Q 
(see Fig. [TJd). However, these real- world networks al- 
ready present short-cuts overlaying the fractal structure, 
and therefore rather than a pure power-law, the scaling 
shows a cut-off at large b, like G' in Fig. [TJd. There- 
fore, these networks are composed by a fractal underlying 
structure, analogous to Go, with some short-cuts gener- 
ating the network G'. The question we want to answer 
here is, what is the a exponent of the short-cuts over- 
laying the fractal network? Since one cannot obtain the 
value of a directly from the data (as it is not possible 
to distinguish a priori between the links of the fractal 
structure and the short-cuts links) we infer its value by 
treating the real-world network as the network G' and 
measuring directly the value of A using the RG flow. 

In Fig. @] we show the results of the RG flow to a sam- 
ple of the WWW 0, a protein homology network [Til ], 
the metabolic network of E. coli fl5| . a yeast protein 
interaction network [l6| . and the co-acting network of 
IMDB that have been found to exhibit fractal topolo- 
gies Q. Table U shows a summary of the exponents for 
real- world networks. For instance we find that while the 
WWW exhibits fractal scaling in Nj, [3| it also presents 
enough short-cuts that its structure belongs to the un- 
stable phase. Thus, the WWW is fractal up to a given 
length scale and then it crosses over to small world behav- 
ior at large scales (iVf, presents an exponential cut-off at 
large b, Fig. [TJd). In fact, for large values of b we observe a 
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FIG. 4: Average degree of the renormalized network versus 
xt = No/Nb = b dB for the studied samples. 



TABLE I: Exponents obtained for the real-world networks. 



Network 


ds from [3] 


A (Fig. IK 


s from Eq. © 


phase 


WWW 


4.1 


0.68 


1.32 


unstable 


Metabolic 


3.5 


0.22 


1.78 


unstable 


PIN 


2.2 


-0.07 


2.07 


stable 


IMDB 


3.88 


-0.71 


2.71 


stable 


PHN 


2.5 


-0.83 


2.83 


stable 



deviation from the power-law, which may be attributed 
to the finite size of these networks, similar to what we 
observe in Fig. [Sfc for the fractal model at different sizes. 
The RG determines the crossover scale such that under 
enough RG steps the WWW finally becomes a complete 
graph. Thus, the renormalization allows to conceptu- 
alize the apparent discordance between small-world and 
fractal properties. Also, due to the its proximity to the 
a = ds point, the WWW is sufficiently randomized to 
give a topology close to optimal information flow. On the 
contrary, the biological networks PHN and PIN, and the 
social network of co-acting belong to the stable phase in- 
dicating that the short-cuts are minimal (the metabolic 
network is unstable but close to the transition point). 
The biological networks display a modular deterministic 
structure shaped by evolution which exhibit pure frac- 
tal character that may be seen as a means of protection, 
preservation and conservation. 

In summary, the RG approach finds the type of short- 
cuts in a given network and determines the location of a 
network in the space of configurations. When the expo- 
nent of the short-cuts is a > 2<ie, the network structure 
belongs to the stable phase, where RG exhibits a fixed 
point consisting of a pure fractal network in the space of 
configurations. On the other hand, when a < 2cIb, the 
network is in the unstable phase where short-cuts become 
dominant, changing dramatically the global distance be- 
tween nodes, and leading to a small-world network at 
large scales. 
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